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ABSTRACT 



The extraordinary SN 1997cy associated with GRB 970514 has been observed 
photometrically and spectroscopically for nearly 2 yr. At the time of discovery 
SN 1997cy was the brightest SN ever observed (Mv< -20.1, Vhd = 19140 km 
s^^, Hq = 65 km s^^Mpc^^). Up to the last available observations (600 days 
after the GRB) the total time-integrated flux was equal to or larger than that 
expected from the complete thermalization of the 7-rays produced by 2.3 M^oi 
^^Co. However, starting already on day 60 the luminosity decline is slower 
than the ^^Co decay rate, indicating that the SN ejecta was interacting with 
circumstellar material (GSM). The interaction appeared to weaken around day 
550. The spectra of SN 1997cy are dominated at all epochs by Ha emission, 
which shows at least 3 components of different widths, as in SN 1988Z. Several 
other lines with different widths are also visible, especially at early epochs. The 
entire light curve of SN 1997cy is reproduced by a model of the interaction of 
the very energetic (£" = 5 x 10^^ ergs) ejecta of a massive star (25 Mq) with the 
GSM, with some contribution from radioactive decays. The GSM could have 
been ejected with a mass-loss rate oi M ^ 4 x 10~"^ Mq yr""*^ as the progenitor 
star evolved from a BSG to a RSG about 10*^ yr before the explosion. The lack 
of oxygen and magnesium lines in the spectra at nebular phases poses a problem 
for models requiring high mass progenitors. The possibility that most of the 
core material of the progenitor has fallen onto a massive black hole so that the 
reverse shock dies at the inner edge of the H/He envelope is discussed. 

Subject headings: supernovae: general — supernovae: individual: SN 1997cy — 
Gamma rays: bursts 
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1. Introduction 

The spatial and temporal coincidence between the GRB980425 and SN 1998bw 
( palama et al. 1998| , [Iwamoto et al. 1998|) raised the issue of the association of (some) 
GRBs with SN explosions. So far the study of the statistical correlation of SNe with BATSE 
and BATSE/Ulysses bursts has given contradictory results. Some authors have proposed 
that all GRBs originate from SNe Ib/c. These are core-collapse SNe whose progenitors have 
lost most or all of the hydrogen envelope ( |Wang fc Wheeler 19981 ) . However, only a weak 
correlation was found between the general list of SNe and GRBs ( |Kippen et al. 1998| ). As 
an alternative it has been suggested that only a small fraction of the GRBs originate from 
asymmetric explosions of rare, highly energetic SN Ib/c ( [Iwamoto et al. 1998| ; [Woosley" 



pj^astmann, fc Schmidt 1999|) . In particular. Bloom et al. (1998) on the basis of a model for 
the radio emission of SN 1998bw have derived the characteristics which a SN must have 
in order to produce a detectable GRB (weak single burst). Again they indicate highly 
energetic Ib/c events as best candidate sources of (some) GRBs. 

Germany et al. (1999) proposed GRB 970514 as another compelling association between 
a highly energetic SN and a GRB. SN 1997cy was discovered on July 16, 1997 ( permany et| 
al. 1997 ) in the compact and faint galaxy n.342 of the cluster Sersic 40/6 = Abell 3266 (z 



0.059; a = 1211 km s ^; [Green et al. 199D ). The relatively narrow Ha emission led to 



the classification of SN 1997cy as a Type Iln. The epoch of the SN explosion is constrained 
by a prediscovery limit on March 12 and is consistent, within the uncertainties, with the 
association with GRB 970514, whose 3.7° error box was centered only 0.88° away. 

The probability of a chance association of the SN with one of the 119 BATSE events 
which occurred in the period between the pre-discovery observation and the discovery 
images is 0.7% ([Germany et al. 1999| ). After SN 1998bw, this is the most compelling case 



for SN/GRB association. GRB 970514 was a single-peak burst detected also above 300 
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keV. At the distance of SN 1997cy, the burst energy would be ~ 4 x 10^^ erg, which makes 
GRB 970514 more energetic than GRB 980425, but orders of magnitude weaker than other 
bursts with measured redshifts ( |Germany et al. 1999| ). 



In this Letter we present and discuss optical observations of SN 1997cy obtained up 
to 2.5 years after discovery, together with a model of the light curve derived with a high 
progenitor mass and explosion energy. 

2. Light Curve 

Imaging photometry of SN 1997cy has been obtained in BVRI with several telescopes 
at ESO-La Silla. Typical seeing was 1.2 arcsec. A sequence of local reference stars has been 
measured on seven photometric nights in order to calibrate non-photometric data. Very 
deep VLT/FORSl imaging (seeing 0.6 arcsec) in V and R about 2.5 yr after discovery did 
not detect the SN. Although there is very good agreement for the stars in common with 
Germany et al. (1999), we note a difference between our SN magnitudes and their estimates. 
This is small (0.15mag) at maximum, but it increases with epoch up to 0.4 mag. We verified 
that measuring the SN magnitude with different techniques, PSF-fitting or subtraction of 
pre-discovery images (kindly provided by B. Schmidt, priv. comm.) - the method used by 
Germany et al. - produces similar results even at late epochs. Artificial star experiments 
show that at all epochs for which we could make the comparison between the data sets our 
typical errors are smaller than 0.1 mag. We believe that the inconsistency is partially due 
to the different passbands ( Permany et al. 19991) , although the modest spectral evolution 
of SN 1997cy should not cause a temporal dependence, and/or to incorrect flux calibration 
of low signal-to-noise data. In particular there is the possibility of contamination by two 
nearby (2 arcsec) bright knots (SE and SW) inside the parent galaxy, which could be severe 
for imaging at small scales or with bad seeing. 
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The light curve of SN 1997cy does not conform to the classical templates of SN II, 
namely Plateau and Linear, but resembles the slow evolution of the Type Iln SN 1988Z 
( P?uratto et al. 1993|) . Assuming Ay = 0.00 for the galactic extinction ( Purstein fc Heiles 



|1982|) we get for SN 1997cy an absolute magnitude at maximum Mv< —20.1. This makes 
SN 1997cy the brightest SN II discovered so far, brighter than the very luminous SN 1979C 
and 1983K (My = -19.9 and -19.3, [Patat et aLT^U^, with Hq = 65 km s'^ Mpc"^). 



Note that Schlegel et al. (1998) give Ay = 0.07. As far as we know, no detections at other 
wavelengths have been reported, except for some not very stringent radio limits at 20 cm 
( permany et al. 1999|) . 

The uvoir bolometric light curve of SN 1997cy obtained using the available BVRI 
photometry and the spectrophotometry is shown in Fig. |l|. Integrating the bolometric 
light curve we estimate that the SN has emitted about 2 x 10^° ergs from the day of 
discovery to our last detection, not including the GRB luminosity. We stress that the uvoir 
determinations are by necessity lower limits of the true bolometric luminosities. 

3. Radioactive decay vs. Interaction 

The late time light curve of most SNe is powered by the radioactive decay of ^^Co to 
^^Fe. If the hard 7-rays from the decay are completely thermalised in the ejecta, the mass 
of ^^Ni, the parent isotope of ^^Co, can be estimated directly from the luminosity on the 
tail of the light curve. Typical ^'^Ni masses for SNe II are ~ 0.08 Mq (e.g. SN 1987A), 
but extreme cases range from 0.002 M© for SN 1997D ( [Turatto et al. 1998|) to 0.3 Mq for 
SN 1992am ( [Schmidt 2000|) . 

Between 60 and 120 days after the outburst, which is taken to be coincident with the 
GRB, the light curve of SN 1997cy matches the decline rate expected if the energy input 
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was the decay of ^^Co. The estimated initial mass of ^^Ni is 2.3Mq. At 4 months the SN 
dechne becomes slower, suggesting that another source of energy is present. Between days 
250 and 550 the light curve again matches the ^^Co decay of QMq, which is clearly dubious. 
Then the light curve suddenly drops, as does the Ha emission. On day 655 the SN could 
not be detected with the ESO 3.6m telescope. The last points of our bolometric light curve 
lie above the radioactive tail of 2.3 Mq of ^^Ni, but are compatible with a value close to 
that proposed by Germany et al. (1999) if we allow for some interaction. However, our 
observations indicate that the late drop mentioned above occurs about 100 days later than 
estimated by Germany et al. 

Other observable signatures of the interaction, besides the slow decline, are complex 
line profiles evolving with time, easily detectable in the prominent Ha line, and powerful 
radio and X-ray emission ( [Aretzaga et al. 1999| ). Ha in SN 1997cy shows an unresolved 
emission and at least two broader components, which become narrower with time. At the 
epoch of the first observation the broadest component has FWHM = 12800 km s~^, and 
its flux dominates over the intermediate (FWHM = 4300 km s^^) one. One year after the 
explosion, the broadest component has almost disappeared and the intermediate component 
(now with FWHM = 2000 km s~^) is the most evident spectral feature in the spectrum. 
SN 1997cy was last detected on April 12, 1999, when only a faint Ha emission was detected 
at the 3-ff level (Fa = 2 x 10^^^ ergs s"^ cm^^; FWHM = 1000 km s"^). 

In Figure 2 we show representative spectra of SN 1997cy at various epochs, and in 
Figure 3 we compare the spectra of SN 1997cy and SN 1988Z, since a certain spectral 
similarity between the two SNe is possibly suggestive. Lines of different widths are present 
in the spectra of SN 1997cy. At early epochs there are very broad bands (FWHM ~ 300 
A) probably blends of several lines, the strongest of which are measured at 4150, 4900, 
5800, 7800 and 9100 A (Call IR triplet). Narrow lines include Ha, B.(3, and [OIII] 5007 
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A. Intermediate width lines (FWHM ~ 2500 - 4300 km s~^) include Hel 5876 A, which 
appears in emission between day 100 and 150 (cf. Fig. 3). 

A multiwavelength study of SN 1988Z showed that the energy radiated from radio to 
X-rays was as high as 10^^ ergs (Aretxaga et al. 1999), suggesting complete reprocessing of 
the mechanical energy of the ejecta in only a few years. The behaviour of the interaction 
between the SN ejecta and a dense Circumstellar Medium in SN 1988Z has been modelled 
by Terlevich et al. (1992) in the framework of the Compact Supernova Remnants (cSNR) 
scenario, and by Chugai & Danziger (1994), who invoked the presence of a rarefied wind 
together with a denser component in the form of clumps or an equatorial wind. The optical 
light curve, the X-ray emission, the Ha emission and the line width evolution of SN 1988Z 
are generally reproduced by these models but the resulting ejecta and CSM masses are 
model dependent. In what follows both the similarities and differences between SN 1988Z 
and SN 1997cy should be kept in mind. The two SNe are similar in apparently having a 
3-component profile of the dominant Ha emission line, and their light curves, though not 
equally well sampled, have the same slow decline for the first two years. SN 1988Z was an 
order-of-magnitude less luminous than SN 1997cy, while the complex of line blending in the 
4000-6000 A region is clearly different, as shown by careful inspection of Fig. 3. 

4. Modelling the light curve 

The similarities between SN 1997cy and SN 1988Z suggest that we investigate an 
interaction model for SN 1997cy as well. Our exploratory model considers the explosion of 
a massive star of M = 25Mq with a parameterized explosion energy E. We assume that 
the collision starts near the stellar radius at a distance ri, where the density of the CSM is 
pi, and adopt for the CSM a power-law density profile p oc r". The parameters E, pi, and 
n, are constrained from comparison with the observations. The numerical code and input 



physics are the same as developed by Suzuki & Nomoto (1995). 

As in the cSNR model the regions excited by the forward and reverse shock emit mostly 
X-rays. The density in the shocked ejecta is so high that the reverse shock is radiative and 
a dense cooling shell is formed (e.g. Suzuki & Nomoto 1995; [Terlevich et al. 1992| ). The 



X-rays are absorbed by the outer layers and core of the ejecta, and re-emitted as UV-optical 
photons. Narrow lines are emitted from the slowly expanding unshocked CSM photoionized 
by the SN UV outburst or by the radiation from the shocks; intermediate width lines come 
from the shock-heated CSM; broad lines come from either the cooler region at the interface 
between ejecta and CSM, or the unshocked ejecta heated by inward propagating X-rays. 

Figure ^ shows the model light curve which best fits the observations. The model 
parameters are: E = 5 x 10^^ erg, pi = 4 x 10~^^ g cm~^ at ri = 2 x 10^*^ cm (which 
corresponds to a mass-loss rate of M = 4 x 10"'^ Mq yr~^ for a wind velocity of 10 km 
s~^), and n = —1.6. The large mass-loss episode giving rise to the dense CSM is supposed 
to occur after the progenitor makes a loop in the HR diagram from BSG to RSG. In this 
model, the mass of the low- velocity CSM is ~ 5 Mq, which implies that the transition from 
BSG to RSG took place about 10^ yr before the SN event. 

The large CSM mass and density and the very large explosion energy are necessary to 
have large shocked masses and thus to reproduce the observed high luminosity. In models 
with low E and high pi the reverse shock speed is too low to produce a sufficiently high 
luminosity. For example, a model with E = 10^^ erg and pi as above yields a value of 
LuvoiR lower than the observed one by a factor of ~ 5. For high E or low pi, the SN ejecta 
expand too fast for the cooling shell to absorb enough X-rays to sustain the luminosity. 
Therefore, E and M are constrained within a factor of ~ 3 of the reported values. 

The shape of the light curve constrains the circumstellar density structure, li n = —2, 
the case of a steady wind, is used, Luvoir decreases too rapidly around day 200. The 
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observed light curve shows a first, shght drop after day ~ 300. This can be reproduced if 
the CSM density is assumed to drop abruptly at the radius reached by the forward shock at 
day 300, so that the collision becomes weaker afterwards. This change of the CSM density 
marks the transition of the progenitor from BSG to RSG, which should have occurred ~ 10^ 
yr before the SN explosion. This is consistent with the simultaneous decrease of the Ha 
luminosity mentioned in Sect. 4. 

After this date, the light curve is powered only by the reverse shock propagating into 
the ejecta. After day ~ 550 the observed light curve shows a second, sharper drop, which 
is reproduced (Figure ^ assuming that after the reverse shock has propagated through ~ 
5 Mq of ejecta it encounters a region of very low density and thus it dies. In other words, 
the model assumes that most of the progenitor's core material fell into a massive (~ IOMq) 
black hole, while only the extended H/He envelope of ~ 5 Mq was ejected. Therefore, both 
forward and reverse shocks propagated inside H/He dominated ejecta. 

However, the fact that the ejecta-CSM model with high explosion energy reproduce the 
light curve does not necessarily mean that the radioactive energy input is small, because 
of the uncertainty in the model fitting. To constrain the contribution from the radioactive 
decay, we calculated the optical light curve of the SN caused by the combination of shock 
heating and radioactive decay. We find that a contribution of up to 0.7 Mq of ^^Ni is 
allowed. However, the decline rate of the radioactive component after day 60 is much faster 
than the half-life of the ^^Co decay, because the explosion is so energetic that only a fraction 
of 7-rays is trapped in the ejecta. This implies that radioactive decay cannot explain the 
slope of the observed light curve even between 60 and 120 days. 
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5. Discussion 

The spectra of SN 1997cy obtained within about 100 days from the GRB show broad 
emission lines which disappear at later times leaving only a blue continuum. Broad emission 
lines originate in the fast expanding ejecta which, in the interaction model, absorb X-rays 
produced in the shocks. This is in principle not very different from the mechanism active in 
the case of radioactive powering, i.e. the deposition of hard radiation (7-rays in that case). 
The width of the broad lines, 13,000 km s^^, is consistent with the outer velocity of the 
expanding ejecta in the interaction model (~ 15,000 km s~^). 

Broad Ca II emission lines are easily identified. However, the Fe II] and Fe III] emission 
lines strong in the nebular spectra of Type la SNe (SNe la) are not obviously identified. 
This could mean that either the ejected ^^Fe (i.e. ^^Ni) mass is small or that the densities 
are so high that forbidden lines are collisionally quenched. In fact the broad feature near 
4500 A might be attributed to the Fell multiplets 37 and 38, while the very broad feature 
near 5000A may be dominated by a blend of the Fell multiplets 42, 48, and 49. 

High densities in the ejecta of SN 1997cy appear on the other hand an almost inevitable 
conclusion: with a characteristic velocity not much higher than that of SNe la nebulae, the 
interaction model for SN 1997cy predicts an ejecta mass of ~ 5 M©, i.e. a factor of ~ 4 
larger. The epoch of the SN 1997cy spectra is also much earlier than it is for SNe la in their 
nebular epoch (100 v. 300 days), thus adding another factor of ~ 3^ to the density. Thus 
we can expect the typical density in the SN 1997cy ejecta to be ~ 100 times that of SNe la 
at 300 days. This implies densities of about 10^ cm~^. Electron densities of course depend 
also on the degree of ionization, but values of at least 10^ cm^'^ can be expected. 

If the forbidden lines are quenched, emission could occur in permitted Fe II lines 
which may be present in the spectrum of SN 1997cy, as suggested above. Nevertheless, the 
excitation of Fell lines in SNe and AGN's involves processes not yet fully understood. The 
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Ha profile of SN 1988Z was explained within the ejecta-CSM interaction scenario by Chugai 
& Danziger (1994), who derived a mass for the ejecta smaller than 1 Mq, with a progenitor 
main sequence mass of 8 to 10 Mq. Modelling the light curve of SN 1988Z as the result of 
interaction may therefore be worthwhile, and it might indicate a larger mass for SN 1988Z. 
However, another property SN 1997cy and SN 1988Z have in common is that neither [01] 
6300 nor Mgl] 4571 have been detected at any stage in their spectra. The critical densities 
of these lines are expected to be a few times 10® and a few times 10^, respectively. The 
envelope density at late phases would not be expected to exceed these values. Since only 
stars with masses less than 8-10 Mq (Nomoto 1984) are expected to synthesize little oxygen, 
one is left with the conundrum of the missing oxygen and magnesium if it is assumed that 
these SNe originated from significantly higher mass stars. 

The model described in the previous section might resolve this conundrum: O and Mg 
synthesised in the core have fallen into the black hole; the ejecta are basically the H/He 
layers and thus contain the original (solar abundance) heavy elements plus some heavy 
elements mixed from the core before fall back. If the SN formed a black hole, the progenitor 
mass was probably larger than ~ 25Mq (e.g. Ergma & van den Heuvel 1998). 

The large energy implied from the modelling of SN 1997cy and indeed its high luminosity 
alone might suggest a causal connection between energetic SN events (hypernovae) and the 
occurrence of GRB's. It appears that hypernovae occur in a variety of SN types. Objects 
such as SNe 1998bw and 1997ef (Iwamoto et al. 2000) are of Type Ic, while SN 1997cy is 
of Type Iln. The recent SN 1999E (Cappellaro et al. 1999), whose spectrum is identical to 
that of SN 1997cy, may be another Type Iln hypernova. Objects such as SN 1988Z also 
deserve renewed study for the reasons discussed in this paper. 

This work has been supported in part by the grant-in-Aid for COE Scientific Research 
(07CE2002) of the Ministry of Education, Science, and Culture in Japan. 
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Fig. 1. — The uyoiy (0.32 to 1 /i) bolometric light curve of SN 1997cy (using the luminosity 
distance, Ab = and M^ q = 4.72). Squares come from the BVRI photometry and triangles 
from the spectrophotometry. The dashed lines indicate the ^^Co decline slope. The rest 
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Fig. 2. — Three spectra describing the evolution of SN 1997cy. The ordinate axis refers to 
the first spectrum (top). Other spectra have been shifted downwards by 3 x 10^^^ ergs s^^ 
cm~^ A~^ and the last multiplied by a factor 3. 
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Fig. 3. — Comparison of the spectra of SNe 1997cy and 1988Z (Turatto et al. 1993) 
at comparable epochs. Both objects show broad hnes in the blue part as well as a 
multicomponent profile of Balmer lines. 
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Fig. 4. — The synthetic hght curve for SN 1997cy obtained with circumstellar interaction 
model. The model parameters are: E = 5 x 10^^ erg, pi = 4 x 10^^^ g cm~^ at ri = 2 x 10^^ 
cm (which corresponds to a mass-loss rate of M = 4 x 10~^ Mq yr~^ for a wind velocity 
of 10 km s^^), and n = —1.6. Shown are the total luminosity from the shocked ejecta Ltot 
(the dotted curve), the UV-optical luminosity -^uvom (solid curve), and the luminosity of 
the X-rays escaping from the ejecta Lx (dashed curve), where Ltot = Luvoir + -^x- 



